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Abstract 

The feasibility of recovering the waste heat from a small-scale incinerator (designed by Industrial Technology Research Institute) and 
generating electric power by a linear free-piston Stirling engine is investigated in this study. A heat-transfer model is used to simulate the 
integration system of the Stirling engine and the incinerator. In this model, the external irreversibility is modeled by the finite temperature 
difference and by the actual heat transfer area, while the internal irreversibility is considered by an internal heat leakage. At a fixed source 
temperature and a fixed sink temperature, the optimal engine performance can be obtained by the method of Lagrange multipliers. 

From the energy and mass balances for the interesting incinerator with the feeding rate at 16 t/d, there is enough otherwise wasted 
energy for powering the Stirling engine and generate more than 50 kW of electricity. 

© 2007 Elsevier Ltd. All rights reserved. 
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1. Introduction 

Incineration with heat recovery is an attractive process 
of recovering energy in some countries [1,2]. The conver¬ 
sion of waste heat to power is very desirable to the areas 
without sufficient energy resources. It has been a common 
practice to convert water into superheated steam by 
absorbing the otherwise wasted energy, the obtained steam 
can then be used to power a steam turbine and generate 
electricity (called steam turbine GENSET). For instance, in 
Japan 769 MW of electric power was generated by 181 
large-scale incineration plants in 1998 [3]. However, it is 
not appropriate to use steam turbine GENSET to recover 
energy from a small-scale incinerator (say less than 5001/ 
d). This is because the thermal efficiency of steam turbine 
will drastically drop as the turbine capacity decreases [4,5]. 
On the other hand, Stirling engine can retain high efficiency 
over power range from a few tens of watts to several 
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kilowatts and is very suitable to use as Stirling GENSET 
for a small or medium scale incinerator. 

Among several types of Stirling machines, the free-piston 
Stirling engine matched with a linear alternator (called 
free-piston Stirling GENSET) has the most mechanical 
simplicity and the least friction loss due to mechanical 
movement [6]. Some technical difficulties inherent to linear 
alternator such as the piston’s overstroking and under¬ 
damped problems have been gradually resolved [7]. With 
many great features in power generation, the free-piston 
Stirling GENSET has been considered in this study as the 
potential candidate to recover the waste energy from a 
small-scale incinerator, designed by the Industrial Tech¬ 
nology and Research Institute (ITRI) of Taiwan. 

Since the Stirling engine is an external combustion 
engine, which means the energy used to power the engine 
movement is entirely obtained from the external source, the 
losses due to the heat transfer in and out of engine or due 
to the heat bypass within the engine structure may be more 
important than that of mechanical losses. Recently, 
Kongtragool and Wongwises [8,9] performed the studies 
and found the power and efficiency of a Strling engine were 
strongly affected by the irreversible heat transfer process 
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Nomenclature 

^source 

heat transfer from the heat source (W) 



Qt 

thermal leakage and other internal irreversibil- 

An,H 

the inner area of engine’s hot-end surface (m ) 


ities (W) 

^oukC 

# ^ 
the outer area of engine’s cold-end surface (m ) 

T c 

the surface temperature of engine’s cold-end 

^out,H 

the outer area of engine’s hot-end surface (m ) 


(K) 

^H,0 

the view factor from hot-end surface to source 

Tfgas 

the temperature of flame gases (K) 


surface 

t h 

the surface temperature of engine’s hot-end (K) 

^coolant 

the convective coefficient of cooling water (W/ 

Tc 

lowest fluid temperature (K) 


m 2 /K) 

^~sink 

the temperature of heat sink (K) 

^He 

the convective coefficient of helium gas (W/m / 

T 

± source 

the temperature of heat source (K) 


K) 

Tc 

highest fluid temperature (K) 

(HA) c 

the thermal conductance from the inner fluid to 

w 

cycle-average power (W) 


the outer cold-end surface (W/K) 

a 

absorptivity of the surface exposed to flame 

(HA) fga8 

s the thermal conductance between hot-end 

<5 

(HA) smk /(HA) c 


surface and flame gases (W/K) 

8 

(HA)c/(HA)h 

(HA) H 

the thermal conductance from the outer hot- 

0 Carnot 

the thermal efficiency of Carnot cycle 


end surface to the inner fluid (W/K) 

Oendo 

the thermal efficiency of internal reversible 

(FT4) sink the thermal conductance from the outer 


cycle 


cold-end surface to the heat sink (W/K) 

f] 

thermal efficiency 

Qc 

the heat transfer from working fluid to engine’s 

9 

(HA) smk /(HA) H 


cold-end surface (W) 

} 

% 

eigen-values of Lagrange Multiplier equations 

Qconv 

convection heat transfer from heat source to 

z 

the ratio of heat which actually flow into the 


hot-end of engine (W) 


fluid 

Qh 

heat transfer from hot-end surface to working 

G 

Stefan-Boltzmann constant, 5.67 x 10 _8 W/ 


fluid (W) 


m 2 K 4 

Srad 

radiation heat transfer from heat source to hot- 

X 

a A H <>F H / (HA) H 


end of engine (W) 

*Pm 

constraint function m 

£?sink 

the rejected heat transfer from working fluid 




(W) 




[8], and also affected by regenerator effectiveness and dead 
volume [9]. 

In this study, a heat transfer model is used to simulate a 
realistic Stirling engine, the effects due to engine’s 
irreversibilities on the output power and thermal efficiency 
will be modeled by mathematical formulation, and 
Lagrange multipliers method is employed to optimize 
the net-work output under the fixed source and sink 
temperatures. 

2. Incinerator and waste heat 

The designed incinerator is of a modular type composed 
of a feeding gear, a dual combustor with primary and 
secondary combustion chambers, a heat recovery system, 
and an air pollution control system. The burning capacity 
ranges 16-50 t/d. 

In order to comply to the regulation set by the 
Environmental Protection Administration of Taiwan to 
eliminate toxic gases, the temperature of secondary 
chamber has to be kept above 1273 K. Therefore, the 
Stirling GENSET is intended to be installed on the 
boundary of secondary combustor with the hot-end surface 
embedded inside the chamber, see Fig. 1. The control of 
chamber’s temperature is automatically made by adjusting 


an auxiliary fuel pump, the flow rate of preheated air from 
the heat recovery system, and the air flow rate directly from 
outside ambient. 

The available waste heat at the feeding of 16 t/d was 
estimated to be 314kW from the energy balance for both 
combustors [10]. The estimation was performed under the 
conditions of burning the sample waste with 100% excess 
air. A complete combustion was assumed, lower heating 
values were used in the calculation, and no auxiliary fuel 
was burnt. The content data for sample waste was provided 
by ITRI, in which they have analyzed the buck wastes 
collecting from four different industrial parks. 

3. Heat transfer model 

To recover the waste heat, the hot-end surface of Stirling 
engine is embedded inside the combustion chamber, see 
Fig. 1. The temperature of heat source T source , shown in 
Fig. 2, represents the controlled chamber temperature of 
incinerator. The temperature of heat sink T sink in Fig. 2 
stands for the temperature of coolant used to cool down 
the cold-end surface, see Fig. 1 . T H is the surface temperature 
of engine’s hot end (receiving flame heat), and T c the 
surface temperature of engine’s cold end (contacting with 
coolant). Since the power cycle (thermodynamic cycle) is 
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completed by the heating, expansion, cooling, contraction 
processes of working fluid, the highest fluid temperature 
(the fluid receiving hot-end heat), 7jj, is lower than T n and 
the lowest fluid temperature (the fluid adjacent to cold-end 
surface), T L , is higher than T c . The magnitude sequence 
of these temperatures is thus to be T source > T H > 
T\J > Ty > T c > T sink . 

3.1. Heat transfer to engine 


temperature of flame gases, we assume 7f gas = T source for 
simplicity purpose. 

3.2. Heat transfer to working fluid 

The cycle-average heat transfer from hot-end surface to 
working fluid is 

Qh = (HA) h (T u - Tv). (4) 


The cycle-average heat transfer from the heat source to 
the engine’s hot end can be expressed as the sum of 
radiation heat transfer and convection heat transfer; i.e. 


^source Srad ~f Sconv’ 

(i) 

where 


Crad = «^out > H^H,0(T(^o u rce - T ti) 

(2) 

and 


2conv = (/T4)f g as (Efg a s — Th). 

(3) 


In Eq. (2), a is the absorptivity of the surface exposed to 
flame and o is the Stefan-Boltzmann constant, which is 
5.67 x 10 8 W/m 2 K 4 . It should be noted that the recipro¬ 
city relation for view factors ( A\F X2 = A 2 F 2 \) has been 
applied in Eq. (2), therefore, ^4 out H is the external area of 
engine’s hot-end and T H ,o is the view factor from hot-end 
surface to source surface. In Eq. (3), (ET4) fgas is the product 
of area ^4 H ,o and the convective heat transfer coefficient 
between hot-end surface and flame gases. T fgas is the 


(ET4) h is the reciprocal of the total thermal resistance from 
the outer hot-end surface to the inner fluid. The thermal 
resistance of a conducting metal (usually Inconel alloy) is 
usually much lesser than the convection resistance of the 
working fluid (usually helium gas). (HA) n can be simplified 
to h He A inH , where /z He is the convective coefficient of 
helium gas and A in n is usually fabricated into a fin shape 
(microfin, etc.) such that A in H is about 3^7 times of v4 out H- 


3.3. Heat rejection from working fluid 

The cycle-average heat transfer from working fluid to 
engine’s cold-end surface is 

Q c = (. HA) C (T L - T C ), (5) 

similar to ( HA)n , (HA) c = hn e A inC where /z He is the 
convective heat transfer coefficient of working fluid 
(helium gas), and A in ^ c is the inner area of engine’s cold- 
end surface. 
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Hot reservoir 



Stirling engine 




Cold reservoir 


Fig. 2. Heat transfer model. 


The rejected heat from working fluid Q c is eventually 
carried away by the cooling water, which is 

2sink = (^^)sink(^L — ^sink)- (6) 

Again, (HA) sink ^A ouUC h coolant where A ouUC is the outer 
area of engine’s cold-end surface and h coo i an t is the 
convective heat transfer coefficient of cooling water. 

3.4. Thermal leakage 

Inside a typical Stirling engine, the temperature of hot- 
end surface (constantly kept above Tjj) is much higher than 
that of cold-end surface (constantly staying below T L ). The 
heat leaking through engine structure (including regenerator 
and cylinder wall) without actually flowing into the working 
fluid to contribute to the power production is inevitable. If 
we express this kind of thermal leakage and other internal 
irreversibilities (such as friction heat etc.) as Q T , see Fig. 2. 
The cycle-average power, W, can be written as 

b 7 = ( Qh — Qt) — (Qc — Qt) 

= (Qh - Qc) 

= (HAMTh - Tv) ~ (HA)c(T L - T C ). (7) 


The thermal efficiency for a heat engine operating between 
heat source at r source and heat sink at T sink , the thermal 
efficiency of Carnot cycle (all processes are reversible) is 

F s j nk 

^Carnot 1 • (^) 

source 

If the reality that temperature difference must exist for the 
heat flow from source to working fluid or from working 
fluid to sink are considered the thermal efficiency should be 
written as 


_ W 

hendo 

2h 



T L 

Tv’ 



where ^ en do represents the thermal efficiency of internal 
reversible (or endor-reversible) cycle, and T L and Tjj are 
the lowest and highest temperatures of engine’s working 
fluid. In reality T source > Tu and T k > r sink , all the heat 
transfer processes are irreversible, and because the energy 
transfers are conducted across the engine’s boundary, these 
processes are called external irreversibilities. 

In addition, when the internal irreversibility is consid¬ 
ered as the internal heat-leakage Q T , the amount of heat 
actually flow into the working fluid (i.e. actually contribute 
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the power generation) is ( Q R — Q T ), we can express 


W = l T l 
2h — Qt 



Here, Q T is used to represent the overall internal 
irreversibility. If we define the ratio of heat which actually 
flow into the fluid as £, i.e., 


£ = 


Qu ~ Qt 

Qh 


the thermal efficiency of a more realistic engine, Eq. (10). 
can be rewritten as 


W Ty 

' , = ¥u = t { '-r 0 . 


(ii) 


4. Optimization (Lagrange multiplier method) 

From designer’s point of view, it is interesting to know 
what kind of temperature combination will result in the 
best performance if the engine is operating under the fixed 
source and sink temperatures. 

In order to obtain the best performance, the work done 
and represented by Eq. (7) is considered as the objective 
function to be maximized. Eq. (7) is rewritten as the 
objective function of 

W = (HA) u (T n -Tv- sT l + sTc), (12) 

where e = (HA) C /(HA) h 

Since the maximum power is limited by Eq. (11) for an 
actual engine, Eq. (11) is rewritten and used as the first 
constrained function: 

(HA) u (T h - TV- sT h + sT c ) = / _Tf\ 
(HAUTu-Tv) + T v j' 

or the constraint function 1 : 

Tc-T l +- s (Th-T v ) 

_t JTu-TQ(T H -T v) = 0 
8 T u 

In addition, the power of engine is closely related to the 
heat flows in and the heat flows out of engine. The second 
constraint is Q u = 2 source , from Eqs. (l)-(4). 

a^H,0^H,0 O' (^source ~ ^h) + (7L4) H s( ^source “ ^h) 

= (HAMTh - Tu), 
or the constraint function 2\ 

T T 4 somce - T T 4 h + 0 T source - 0 T H - T H + T v = 0, (14) 

where z = aA^oF^oofHA)^ and 6 = (HA) sink /(HA) u . 
The third constraint is Q c = 2 sink , from Eqs. (4) and (5). 

(HA) c (T l - Tc) = ( HA) sink (T c - r smk ), 


or the constraint function 3: 

T l = T c + §T c - ST sink , (15) 

where 5 = (HA) sink /(HA) c . 

If Source and T sink are fixed, and the thermal con¬ 
ductances are set to be some practical values, the best 
combination of T n , Tu, T L , and T c can be obtained by 
solving the Lagrange multiplier equations, which are: 

3 

VW-^2l m V(f> m = 0, (16) 

m =1 


T v , T l , T c ) = T c - T h +fir H - 7u) 

--^Tu-Tu-Afj+Tjj =0 , 

( 17 ) 

02<Th, Tu, T l , T c ) = ource - + 6T somce 

- 6T h - T h + Tu = 0, (18) 

Tu, T h , T c ) = T l -T c - ST c + ST smk = 0. (19) 

The above constraint functions </>i, <j 2 , and </> 3 , are 
essentially the same as Eqs. (13)—(15). Since there are four 
independent variables (T n , Tu, T k , T c ), the gradient W 
and gradient <j) m shown in Eq. (16) are the derivative of 
objective function and constraint functions to each 
independent variable, which will result in four equations. 
For example if one takes the partial derivative of Eq. (16) 
with respect to Tu, 0(16)/0T H , he gets 

(™) H - M‘-«+«£) 

+ 4zj 2 Tu T + 2-2 = 0. (20) 

Similarly, from 0(16)/07jj, 0(16)/0T L , and 0(16)/0T C one 
will, respectively, obtain the other three equations: 


~(HA) U + 2 



£ , Z T k tA 

s s T 2 V ) 




s(HA) u + 2i (l + - — --=—\ - 23—0 

V e zT\jJ 



8(HA)u - 2i + 2 3 + 2 3 (5 = 0. (23) 

Equations (17)—(23) can be solved simultaneously to obtain 
the optimized value of T H , Tu, T L , T c and the eigen-values 
2 i- 2 3 . The power output, thermal efficiency and other heat 
transfer quantities can then be calculated. 


5. Results and discussion 


It is interesting to know how much of power output can 
be obtained from a properly designed Stirling engine with 
an unit outer hot-end surface, i.e. we set ^4 0 ut,H = 1 m . 
Since (HA)u = hu e A in u, the convective heat transfer 
coefficient of helium gas is about twice of flue gases, and 
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the inner finned surfaces at both hot and cold ends are 
assumed five times of the outer surfaces. It is reasonable 
to set {HA) H = 1000 W/K, {HA) fgas = 100W/K, T source = 
1273 K, and F sink = 300 K. 

For an endor-reversible heat engine operating steadily 
between two fixed temperatures, Bejan [11] found the 
engine’s best performance can be obtained if the total 
thermal conductance ( HA ) is evenly divided for hot-end 
and cold-end heat exchangers, i.e. {HA)u = {HA) C . In our 
previous study [12], we derived a more general {HA) 
formulation for a realistic engine. Under a fixed total 
thermal conductance constraint, (7T4) H is smaller but very 
close to {HA) C if the internal irreversibility is considered. 

From the above reasonings, we set s = {HA) c /{HA)u = 1, 
and d = {HA) s [ n y/{HA) c = 1. When Qj/Qh = 0.3 (or 
cj = 0.7), the solutions of Eqs. (17)—(23) are T u = 995 K, 
T v = 887.9 K, F l = 438.3 K, T c = 369.2 K, = 1203, 
2 2 = —172.1, and 2 3 = 101.6. The resulting heat transfers 
Of Gsomce = 107 > 100W, Qt = 32,130 IF, g sink = 69,200 W, 
network IF = 37,969 W, and thermal efficiency r\ — 0.3544. 
Fig. 3 shows the calculated Q som ce , IF and r\ when the degree 
of heat bypass is changed, and Fig. 4 depicts the associated 
temperature distributions when heat bypass is varied. 
At £ = 1.0 (no heat leakage), the calculated thermal 
efficiency = 53.33%, which is close to the Curzon-Ahlborn 
efficiency of rj CA = l-(T sink /T source f 5 = 51.45%. As ex- 
pected, the thermal efficiency as well as the output power 
decrease when the degree of leakage increases (£ decreases). 
At £ = 0.75 (25% heat bypass), the calculated thermal 
efficiency is about 38.2% which is close to 50% of the Carnot 
efficiency = l—(r sink /r source ) = 76.43%. Since the optimal 
temperatures are slightly increased when the degree of heat 
bypass increases (£ decreases), see Fig. 5, the change rate of 
efficiency is not exactly proportional to the change rate of 
heat bypass. 

If the source and sink temperature are still kept at 1273 
and 300 K, but the value of thermal conductance, {HA ) n , is 
changed, and the ratios of thermal conduction are kept 



Fig. 3. Heat input, network and efficiency at different degrees of heat loss. 



Fig. 4. F h , Tu, T l , and T c values at different degrees of heat loss. 



(HA) h (W/K) 


60% 


40% 


20 % 


sr 


Fig. 5. Heat input, network and efficiency at different 


unchanged, the calculated Q source , IF and rj are presented in 
Fig. 5 and the temperature distributions are presented in 
Fig. 6. Results indicate 2 source , IF and r\ are increasing and 
all of the temperatures are decreasing as the thermal 
conductance is getting larger. Since T u represents the 
highest temperature of engine body, the predicted value of 
T n provides important data to select the hot-end material. 
Fig. 6 indicates the curve of T H has the largest decreasing 
rate. Since temperature is the easiest data to measure, from 
the designer’s point of view, all of the calculated optimal 
temperatures can be used as reference values to compare 
with the measured data during the development stage. This 
kind of comparisons can provide some insight information 
to perform structural adjustments or local fine turnings to 
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improve the overall design. From Fig. 5, we can see even at 
(HA)n = 500 W/K and £ = 0.7(30% bypass heat), a very 
conservative estimation, the net-work (power output) is 
about 26 kW, and the required input heat is about 75 kW. 
Consider the physical dimensions of secondary combustor, 
there is enough space to install 2 units of Stirling GENSET, 
and there are more than enough waste heat to provide the 
engine’s input heat even at the feeding rate of 16 t/d. 

6. Conclusion 

A simple heat transfer model is used to model the 
inherited irreversibilities of a realistic Stirling engine. The 
heat transfers from heat source to engine’s hot end, and 
from engine’s cold end to heat sink are used to represent 
the external irreversibilities, and the degree of heat bypass 
is used to simulate the overall internal irreversibility. The 
calculated thermal efficiency without internal bypass heat 
compares well with the Curzon-Ahlborn efficiency. At 
25% bypass heat, the calculated thermal efficiency is very 
close to one-half of the Carnot efficiency. 

With 1 m of hot-end surface embedded inside the 
secondary combustion chamber, a well designed Stirling 
GENEST should be able to generate more than 25 kW of 


power. Considering the configuration constraints of the 
secondary chamber and the stability of operating tempera¬ 
ture, more than 50 kW of power can be recovered from the 
interesting incinerator. 
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